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Abstract 

Here, we investigate the growth of matter density perturbations as well as the generalized second 
law (GSL) of thermodynamics in the framework of f{R)- gravity. We consider a spatially flat FRW 
universe filled with the pressureless matter and radiation which is enclosed by the dynamical 
apparent horizon with the Hawking temperature. For some viable f(R) models containing the 
Starobinsky, Hu-Sawicki, Exponential, Tsujikawa and AB models, we first explore numerically 
the evolution of some cosmological parameters like the Hubble parameter, the Ricci scalar, the 
deceleration parameter, the density parameters and the equation of state parameters. Then, we 
examine the validity of GSL and obtain the growth factor of structure formation. We find that for 
the aforementioned models, the GSL is satisfied from the early times to the present epoch. But in 
the farther future, the GSL for the all models is violated. Our numerical results also show that for 
the all models, the growth factor for larger structures like the ACDM model fit the data very well. 
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I. INTRODUCTION 


The observed accelerated expansion of the universe, as evidenced by a host of cosmolog¬ 
ical data such as supernovae la (SNela) [1], cosmic microwave background (CMB) [2, 3], 
large scale structure (LSS) [4], etc., came as a great surprise to cosmologists. The present 
accelerated phase of the universe expansion reveals new physics missing from our universe’s 
picture, and it constitutes the fundamental key to understand the fate of the universe. 

There are two representative approaches to explain the current acceleration of the uni¬ 
verse. One is to introduce “dark energy” (DE) [5] in the framework of general relativity 
(GR). The other is to consider a theory of modified gravity (MG), such as f(R ) gravity, 
in which the Einstein-Hilbert action in GR is generalized from the Ricci scalar R to an 
arbitrary function of the Ricci scalar [6]. Here, we will focus on the later approach. 

In [7], it was shown that a f(R) model with negative and positive powers of Ricci curvature 
scalar R can naturally combine the inflation at early times and the cosmic acceleration 
at late times. It is actually possible for viable f(R) models for late time acceleration to 
include inflation by adding R 2 term. Therefore, it is natural to consider combined f(R) 
models which describe both primordial and present DE using one f(R) function, albeit one 
containing two greatly different characteristic energy scales [8, 9]. In [10], it was pointed 
out that the f(R )-gravity can also serve as dark matter (DM). In [11], a set of f(R )-gravity 
models corresponding to different DE models were reconstructed. Although a great variety 
of f(R) models have been proposed in the literature, most of them is not perfect enough. 
An interesting feature of the f(R) theories is the fact that the gravitational constant in 
f(R)- gravity, varies with length scale as well as with time [12]-[18]. Thus the evolution of 
the matter density perturbation, S m = 5p m /p m , in this theory is affected by the effective 
Newton coupling constant, G e fj, and it is scale dependent too. Therefore, the matter density 
perturbation is a crucial tool to distinguish MG from DE model in GR, in particular the 
standard ACDM model. 

On the other hand, the connection between gravity and thermodynamics is one of sur¬ 
prising features of gravity which was first reinforced by Jacobson [19], who associated the 
Einstein field equations with the Clausius relation in the context of black hole thermody¬ 
namics. This idea was also extended to the cosmological context and it was shown that the 
Friedmann equations in the Einstein gravity [20] can be written in the form of the first law 
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of thermodynamics (the Clausius relation). The equivalence between the first law of thermo¬ 
dynamics and the Friedmann equation was also found for /(A)-gravity [21]. Besides the first 
law, the generalized second law (GSL) of gravitational thermodynamics, which states that 
entropy of the fluid inside the horizon plus the geometric entropy do not decrease with time, 
was also investigated in /(A)-gravity [22], The GSL of thermodynamics in the accelerating 
universe driven by DE or MG has been also studied extensively in the literature [23]-[33]. 

All mentioned in above motivate us to investigate the growth of matter density per¬ 
turbations in a class of metric f(R) models and see scale dependence of growth factor. 
Additionally, we are interested in examining the validity of GSL in some viable /(A)-gravity 
models. The structure of this paper is as follows. In Sec. II, within the framework of /(A)- 
gravity we consider a spatially flat Friedmann-Robertson-Walker (FRW) universe filled with 
the pressureless matter and radiation. In Sec. Ill, we study the growth rate of matter 
density perturbations in /(A)-gravity. In Sec. IV, the GSL of thermodynamics on the 
dynamical apparent horizon with the Hawking temperature is explained. In Sec. V, the 
cosmological evolution of /(A) models is illustrated. In Sec. VI, the viability conditions for 
/(A) models are discussed. In addition, some viable /(A) models containing the Starobin- 
sky, Hu-Sawicki, Exponential, Tsujikawa and AB models are introduced. In Sec. VII, we 
give numerical results obtained for the evolution of some cosmological parameters, the GSL 
and the growth of structure formation in the aforementioned /(A) models. Section VIII is 
devoted to conclusions. 


II. /(A)-GRAVITY FRAMEWORK 


Within the framework of /(A)-gravity, the modified Einstein-Hilbert action in the Jordan 
frame is given by [6] 


Sj 


V~-g d 4 a; 


7(g) 

167tG 


+ L 


matter ? 


(i) 


where G, g, A and L matter are the gravitational constant, the determinant of the metric g jiv , 
the Ricci scalar and the lagrangian density of the matter inside the universe, respectively. 
Also /(A) is an arbitrary function of the Ricci scalar. 

Varying the action (1) with respect to g lw yields 


FG, W = 8t tGT$> - \g^{RF - /) + - g^DF. 


( 2 ) 
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Here F = df/dR, G jiv = R llv — \Rg^ v and is the energy-momentum tensor of the 
matter. The gravitational held equations (2) can be rewritten in the standard form as 


[34, 35] 


G^ = 8nG(T^+T^), (3) 

with 

8t tGTM = (1 - F)G, U - \g^{RF - /) + W,W U F - g^UF. (4) 

For a spatially hat FRW metric, taking T= cliag {—p,p,p,p) in the prefect fluid 
form, then the set of held equations (3) reduce to the modified Friedmann equations in the 
framework of f(R)~ gravity as [36] 

3 H 2 = 8nG(p + p D ), (5) 

2 H = -8nG(p +p D +p + p D ), (6) 


where 


with 


8nGp-D 

8nGpB 


RF-f ) — 3HF + 3H 2 (l — F), 

^-(RF-f) +F + 2HF- (1-F)(2H+ 3H 2 ) , 


( 7 ) 

( 8 ) 


R = 6(H + 2H 2 ). 


(9) 


Here H = ci/a is the Hubble parameter. Also pn and pn are the curvature contribution to 
the energy density and pressure which can play the role of DE. Also p = pbm + Pdm + Prad 
and p = p ra d = Prad/3 are the energy density and pressure of the matter inside the universe, 
consist of the pressureless baryonic and dark matters as well as the radiation. On the whole 
of the paper, the dot and the subscript R denote the derivatives with respect to the cosmic 
time t and the Ricci scalar R. respectively. 

The energy conservation laws are still given by 


Pm + 3 Hp m — 0, (10) 

p rad + 4 F/p rad = 0, (11) 

p D + 3 //(pd + Pd) = 0, (12) 


where p m = pbm + Pdm- From Eqs. (10) and (11) one can find 


P 


Pm 0 Pradp 

O ■ A 


(13) 
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where p mo = pbm 0 + Pdm 0 and p ra d 0 are the present values of the energy densities of matter 
and radiation. We also choose ao = 1 for the recent value of the scale factor. 


Using the usual definitions of the density parameters 
^ _ Pm _ 8nGp mo pj _ Prad _ S^f^Prado 

m = Vc = !uPa?~ , rad = ^T = 3 FT 2 a 4 ’ 


UraH 


Un = ^ = 


p D _ 8irGp D 
p r ~ 3 H 2 ’ 


in which p c = 3 H 2 / ( 87 tG) is the critical energy density, the modified Friedmann equation 
(5) takes the form 

1 — 1 n f^rad T Ud. (15) 

From the energy conservation (12), the equation of state (EoS) parameter due to the cur¬ 
vature contribution is defined as 

Pd . Pd n 

u D = - = -l- ——. (16) 

Pd 3Hpn 

Using the modihed Friedmann equations (5) and ( 6 ), the effective EoS parameter is obtained 


P + Pd , 2 H 

u e & = — -= -1 - 17 

P + Pd 3 H 2 

Also the two important observational cosmographic parameters called the deceleration q and 

the jerk j parameters, respectively related to a and a , are given by [37] 

a _ H „ R . . 

q = ~^lP~~ ~lP~ _ 6 FP’ ^ ^ 

= + ^P = 2 + q + ^p- < 19 > 


III. GROWTH RATE OF MATTER DENSITY PERTURBATIONS 

Here, we study the evolution of the matter density contrast <5 m = 5p m /p m in /(i?)-gravity. 
To this aim, we consider the linear scalar perturbations around a flat FRW background in 
the Newtonian (longitudinal) gauge as 

ds 2 = -(1 + 2 T)d t 2 + a 2 (t)( 1 + 2 <F)da; 2 , ( 20 ) 

with two scalar potentials T and <f> describing the perturbations in the metric. In this 
gauge, the matter density perturbation S m and the perturbation of 6F(R) obey the following 
equations in the Fourier space [38, 39] 

L+(2H+^\s m -^^S m =-^ (-6H 2 + ^pj5F + 3HSF + 35F , ( 21 ) 
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5F + 3 HSF +(~F^r-^]sF = ^ p m S m + FS m , (22) 

\ a 1 3F r 3 / 3 

where k is the comoving wave number. For the modes deep inside the Hubble radius (i.e. 
k 2 /a 2 3> H' 2 ), we have |F| <C HF and SF <C HSF <C H 2 , hence the evolution of matter 
density contrast S m reads [40, 41] 

5 m + 2 H5 m — AnG e sp m Sm — 0, (23) 


where 

4 _ 1 MV 1 

3 3 V + M 2 a 2 \ 1 

and M 2 = The fraction of effective gravitational constant to the Newtonian one, i.e. 

Gefr/G, is defined as screened mass function in the literature [16]. Equation (24) obviously 
shows that the screened mass function is the time and scale dependent parameter. 

With the help of new variable namely g(a) = S m /a which parameterizes the growth of 
structure in the matter, Eq. (23) becomes 


G fi ff — 


G 

F 


d 2 g 

d In a 2 


+ 



d g 

din a 


+ 



4vrG e frp n 
H 2 


9 = 0. 


(25) 


In general, there is no analytical solution to this equation. But in [42] for an asymptotic 
form of viable f(R ) models at high curvature regime given by f(R) = R + R~ n where 
n > —1, an analytic solution for density perturbations in the matter component during the 
matter dominated stage was obtained in terms of hypergeometric functions. In what follows 
we solve the differential equation (25), numerically. To this aim, the natural choice for the 
initial conditions are g(a m ) = 1 and | a =a m = 0, where a m = 1/(1 + z m ) should be taken 
during the matter era, because for the matter dominated universe, i.e. H 2 = 8nGp m /3 and 
Gefr/G = 1, the solution of Eq. (23) yields 5 m = a. The growth factor is defined as [43] 


/(*) 


d In <5 m 
d In a 


«(! + *) 


d In <5 m 
d^ 


(26) 


which is an observational parameter. In the present work, we obtain the evolution of linear 
perturbations relevant to the matter spectrum for the scales; k = 0.1,0.01,0.001 h Mpc -1 , 
where h corresponds to the Hubble parameter today. 
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IV. GENERALIZED SECOND LAW OF THERMODYNAMICS 


Here, we are interested in examining the validity of the GSL of gravitational thermody¬ 
namics for a given /(A) model. According to the GSL, entropy of the matter inside the 
horizon beside the entropy associated with the surface of horizon should not decrease during 
the time [20]. As demonstrated by Bekenstein, this law is satisfied by black holes in contact 
with their radiation [44], The entropy of the matter containing the pressureless matter and 
radiation inside the horizon is given by the Gibbs’ equation [23] 

T A dS = d E + pdV. (27) 


Taking time derivative of Eq. (27) and using the energy equations (10)-(11) as well as the 
Friedmann equations (5)-(6) one can find 

TA S=^(f A -Hr A )(-2H + H±-Q F , (28) 

where r A = ( H 2 + ^-) _1 / 2 and T A = ^^-(l — -mr A ) are the dynamical apparent horizon 
and Hawking temperature, respectively. The horizon entropy in /(A)-gravity is given by 
‘-A = IF [45], where A = is the area of the apparent horizon. Taking the time 

derivative of S' A one can get the evolution of horizon entropy as 

r A s A = 5 V( 2 ^ A - rA )(^ + 4) f . (29) 

Now we can calculate the GSL due to different contributions of the matter and horizon. 
Adding Eqs. (28) and (29), one can get the GSL in /(A)-gravity as [22] 


T A Stot = 


4 GH 4 L 


2H 2 F - HHF + 2 (H + H 2 )F 


(30) 


where S tot = S + Sa■ Note that Eq. (30) shows that the validity of the GSL, i.e. > 0, 

depends on the /(A)-gravity model. For the Einstein gravity (A = 1), one can immediately 
find that the GSL (30) reduces to 


T A S t ot = 


Ft' 2 


2GH 4 


> 0 , 


(31) 


which shows that the GSL is always fulfilled throughout history of the universe. 
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V. COSMOLOGICAL EVOLUTION 


Here, we recast the differential equations governing the evolution of the universe in di¬ 
mensionless form which is more suitable for numerical integration. To do so, following [46] 
we use the dimensionless quantities 

i=H 0 t, H=^~, (32) 

-Ho -H 0 

f f i ? Z7 T? jji RrR /'qq\ 

J ~ 771’ ^R-77=2> ™ ~ 77^4’ (3d) 

-tt 0 -H 0 ui 0 

where 77o is the Hubble parameter today. With the help of the above definitions and using 

S = -*0+4- < 34 > 

one can rewrite the modified Friedmann equation (5) as follows 

H 2 = n mo [(l+z) 3 +x(l+z) 4 ]+(F-l)[H 2 -(l+z)HH']~(f-R) + (l+z)H 2 F R R ', (35) 

where x — Prad 0 /Pm 0 = n ra d 0 /^m 0 and prime denotes a derivative with respect to the 
cosmological redshift z — - — 1. 

To solve Eq. (35) we introduce new variables as [47]: 

m-= — = ■§-—(i +-) 3 - x(i +-) 4 , (36) 

P mo ^ ^mo 

and 

Vr-= W- 3 (1 + z) 3 . (37) 

^ ^mo 

Taking the derivative of both sides of Eqs. (36) and (37) with respect to redshift z yield 

-(! + z )y'u = ^VR - %h, (38) 

-(1 + z)y' R = 9(1 + z) 3 - W { 9 „ + H-(/ - R) 

-(F-l) ^ - j/h - 4 ((1 + ^) 3 + 2x(l + ^) 4 ) }. (39) 

Finally, inserting Eq. (39) into the derivative of Eq. (38) gives a second differential equation 
governing yn(z) as [48] 


(1 + z) 2 y'ft + Ji(l + z)y' n + J 22 /H + ^3 — 0, 


(40) 



where 


Ji = -3 - 


(1-F 


Jo = 


\6H 2 F R 
2 -F 


3H 2 F r ’ 


J 3 — —3(i + zy 


6H 2 F } 


R 


(1-F) (l + ^) 3 + 2 X (l + ^) 4 + 


311 


<R-f) 


m 0 


(41) 

(42) 

(43) 


Equation (40) cannot be solved analytically. Hence, we need to solve it numerically. To do 
so, we use the two initial conditions yn{z\) = 3 and y' n (z\) = 0 which come from the ACDM 
approximation of f(R) model in high curvature regime. Notice z\ is the proper redshift in 
which we have RF&(zi) < 1 CT 13 . 

With the help of Eqs. (14), (16), (17) and (36) one can obtain the evolutionary behaviors 
of the matter density parameter, H m (z), DE density parameter, Ho ( 2 ), EoS parameter of 
DE, uj-d{z)i and effective EoS parameter, oj e s(z), in terms of yu and its derivatives as follows 

(1 + ^) 3 




[z = 


Hd(^) — 


y H + (l + -) 3 + y(l + -) 4 ’ 

_ Un _ 

Vn + (1 + ^) 3 + x(l + zy 1 

1 + z 


ujd(z) = -1 + 


Vn 

yn 


u e s{z) — —1 + 


(1 + z) 


y'n + 3(1 + z) 2 + 4y(l + z) c 


(44) 

(45) 

(46) 

(47) 


3 L m + 0- + -) 3 + x(i + zy 

Also from Eqs. (18), (19) and (36) one can get the evolutions of the deceleration and jerk 
parameters as 


q(z) = -1 + 
j(z) = 1 + 


(l + z) |~ 2/h + 3(1 + z) 2 + 4y(l + z) 2, 
2 [ yn + (1 + zf + y(l + z) 4 

(1 + *) f (1 + z)Vk ~ 2 y'n + 4 *(1 + z f 

2 L yn + (1 + z f + x(i + zy 


(48) 

(49) 


VI. VIABLE f(R)-G RAVITY MODELS 

Since we are intersected in investigating the growth of structure formation and examin¬ 
ing the GSL in /(F)-gravity, hence in what follows we consider some viable f(R) models 
including the Starobinsky, Hu-Sawicki, Exponential, Tsujikawa and AB models. 
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A. Starobinsky Model 


The Starobinsky f(R) model is as follows [40] 

f(R) = R + XR 


, R 2 

+ m 


(50) 


where n > 0, A and R s are constant parameters of the model. Following [53], we take n — 2 
and A = 1. Note that in the high z regime (z ~ Zi) we have R/R s S> 1. This yields the 
f(R) model (50) to behave like the ACDM model, i.e. f(R) = R — 2A. Consequently, the 
constant parameter R s is obtained as R s = 18f2 mo i/o/A. 


B. Hu-Sawicki Model 


This model was reconstructed based on the local observational data and presented by Hu 
and Sawicki [47] as 


W) 


<>*(£)■ 
<*(£)’ + 1 ’ 


(51) 


where n > 0, ci,C 2 and R s are constants of the model. For this model we take n = 4, 
ci = 1.25 x 10 -3 , c 2 = 6.56 x 10” 5 [46], and obtain R s = 18c2fl mo -^o/ c i- 


C. Exponential Model 

This model is defined by the following function [48], 

f{R)=R^J3R a ( 1-e"^), (52) 

where (3 and R s are two constants of the model. Here R s corresponds to the characteristic 
curvature modification scale. Here we take f3 = 1.8 [48] and obtain R s = 18f) mo i?o//3. 


D. Tsujikawa Model 

This model was originally presented in [39] as 

f{R) = R-\R B tanh(Jpj, (53) 

where A and R s are the model parameters. For this model we obtain R s = 18fl mo f7g/A and 
set A = 1 [54], 
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E. AB Model 


This model was proposed by Appleby and Battye [8, 55] as 


m 


R e 

-h - log 

2 2 & 


cosh — 6) 
cosh(6) 


(54) 


where b is a dimensionless constant and e = R s / [b + log(2 cosh 6)]. The constant R s can be 
obtained at high curvature regime when the AB f(R) model (54) behaves like the ACDM 
model, i.e. f(R) = R — 2A. This gives 


iL 


Here, we also set b = 1.4. 


—36 O mo // ( j b + log(2 cosh b) 
log fbf^j 


VII. NUMERICAL RESULTS 

Here to solve Eq. (40) numerically, we choose the cosmological parameters H mo = 0.24, 
ff Do = 0.76 and fl rado = 4.1 x 10~ 5 . As we have already mentioned, we use the two suitable 
initial conditions yn(z i) = 3 and y'n(zi) = 0, in which z\ is obtained where RFr —> 10~ 13 . 
For the Starobinsky, Hu-Sawicki, Exponential, Tsujikawa and AB f(R) models we obtain 
Zi =15.61, 13.12, 3.66, 3.52 and 3.00, respectively. 

In addition, to study the growth rate of matter density perturbations, we numerically 
solve Eq. (25) with the initial conditions g(z m ) = 1 and (d^/dlna)| Zm = 0, in which z m is 
obtained where fl m (z m ) = 1. For the aforementioned models we obtain z m =14, 13, 12, 14 
and 14.36, respectively. 

With the help of numerical results obtained for yu{z) in Eq. (40), we can obtain the 
evolutionary behaviors of H, R , q, £2 m , fl D , Weff, wd arid GSL for our selected f(R) models. 
The results for the Starobinsky, Hu-Sawicki, Exponential, Tsujikawa and AB f(R) models 
are displayed in Figs. 1-5. Figures show that: (i) the Hubble parameter and the Ricci 
scalar decrease during history of the universe, (ii) The deceleration parameter q varies from 
an early matter-dominant epoch (q = 0.5) to the de Sitter era (q = —1) in the future, as 
expected. It also shows a transition from a cosmic deceleration q > 0 to the acceleration 
q < 0 in the near past. The current values of the deceleration parameter for the Starobinsky, 
Hu-Sawicki, Exponential, Tsujikawa and AB f(R) models are obtained as q 0 = —0.56, 
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—0.60, —0.56, —0.57 and —0.60, respectively. These are in good agreement with the recent 
observational constraint go = — 1 0.43^o;u (68% CL) obtained by the cosmography [56]. (iii) 
The density parameters Dd and Q m increases and decreases, respectively, as z decreases, (iv) 
The effective EoS parameter, o; e g, for the all models, starts from an early matter-dominated 
regime (i.e. u e R = 0) and in the late time, z —>• —1, it behaves like the ACDM model, 
ta e ff —>• — 1. (v) The EoS parameter of DE, cud, for the all models starts at the phase of 
a cosmological constant, i.e. cud = —1 , and evolves from the phantom phase, cud < — 1 , 
to the non-phantom (quintessence) phase, wd > — 1. The crossing of the phantom divide 
line cud = — 1 occurs in the near past as well as farther future. At late times (z —> —1 ), 
ud approaches again to —1 like the ACDM model. Moreover, the present values of wd for 
the Starobinsky, Hu-Sawicki, Exponential, Tsujikawa and AB f(R) models are obtained as 
o;d 0 = —0.94, —0.98, —0.93, —0.94 and —0.97, respectively. These values satisfy the present 
observational constraints [2, 3]. 

(vi) The variation of the GSL shows that it holds for the aforementioned models from 
early times to the present epoch. But in the farther future, the GSL for the Starobinsky, Hu- 
Sawicki, Exponential, Tsujikawa and AB f(R) models is violated for —0.996 < z < —0.955, 
-0.935 < z < -0.909, -0.897 < z < -0.751, -0.997 < z < -0.958 and -0.995 < z < 
—0.950, respectively. To investigate this problem in ample detail, using Eq. (17) we rewrite 
Eq. (30) in terms of aj e g as 


T A S tot 


1 

AG 


— (1 + ui e ^yF + — (1 + ca e ff) 


F 

H 


(1 + 3ta e ff) 


F 
H 2 


(55) 


which shows that in the farther future z —> — 1 when ta e ff —> —1 (see Figs. 1-5) we have 

F 


T A S tot ~ 


2 GH 2 


(56) 


According to Eq. (56), the validity of GSL, i.e. T A S tot > 0, depends on the sign of F. In 
Figs. 1-5, we plot the variation of F/(2F[ 2 ) versus z in the farther future for the selected 
f(R) models. Figures confirm that when the sign of F changes from positive to negative 
due to the dominance of DE over non-relativistic matter then the GSL is violated. Although 
the parameters used for each model in Figs. 1-5 are the viable ones, by more fine tuning 
the model parameters the GSL can be held. For instance, in AB f(R) model by choosing 
the model parameter as b = 1.3, the GSL is always satisfied from early times to the late 
cosmological history of the universe. 
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In Figs. 6-10, we plot the evolutions of RFr, G e s/G, g and the growth factor / versus z 
for the selected f(R) models. Figures show that: (i) RFr goes to zero for higher values of z 
which means that the f(R) models at high z regime behave like the ACDM model, (ii) The 
screened mass function G e g/G for a given wavenumber k is larger than one which makes 
a faster growth of the structures compared to the GR. However, for the higher redshifts, 
the screened mass function approaches to unity in which the GR structure formation is 
recovered. Note that the deviation of G e s/G from unity for small scale structures (larger k) 
is greater than large scale structures (smaller k). (iii) The linear density contrast relative to 
its value in a pure matter model g = 5/a starts from an early matter-dominated phase, i.e. 
g — 1 and decreases during history of the universe. For a given z, g in the all f(R) models, 
is greater than that in the ACDM model, (iv) The evolution of the growth factor f(z) for 
f(R) models and ACDM model together with the 11 observational data of the growth factor 
listed in Table I show that for smaller structures (larger k), the all f(R) models deviate 
from the observational data. But for larger structures (smaller k), the growth factor in the 
all f(R) models, very similar to the ACDM model, fits the data very well. 

VIII. CONCLUSIONS 

Here, we investigated the evolution of both matter density fluctuations and GSL in some 
viable f(R) models containing the Starobinsky, Hu-Sawicki, Exponential, Tsujikawa and AB 
models. For the aforementioned models, we first obtained the evolutionary behaviors of the 
Hubble parameter, the Ricci scalar, the deceleration parameter, the matter and DE density 
parameters, the EoS parameters and the GSL. Then, we explored the growth of structure 
formation in the selected f(R) models. Our results show the following. 

(i) All of the selected f(R) models can give rise to a late time accelerated expansion phase 
of the universe. The deceleration parameter for the all models shows a cosmic deceleration 
q > 0 to acceleration q < 0 transition. The present value of the deceleration parameter takes 
place in the observational range. Also at late times (z —> —1), it approaches a de Sitter 
regime (i.e. q —> — 1), as expected. 

(ii) The effective EoS parameter cu e g for the all models starts from the matter dominated 
era, u; e fr — 0, and in the late time, z —1, it behaves like the ACDM model, toyr —> —1. 

(iii) The evolution of the EoS parameter of DE, cud, shows that the crossing of the 
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phantom divide line ud — — 1 appears in the near past as well as farther future. This is 
a common physical phenomena to the existing viable f(R ) models and thus it is one of 
the peculiar properties of f(R) gravity models characterizing the deviation from the ACDM 
model [54], 

(iv) The GSL is respected from the early times to the present epoch. But in the farther 
future, the GSL for the all models is violated in some ranges of redshift. The physical reason 
why the GSL does not hold in the farther future is that the sign of F changes from positive 
to negative due to the dominance of DE over non-relativistic matter. 

(v) For the all models, the screened mass function G e g/G is larger than 1 and in high 
z regime goes to 1. The deviation of G e g/G from unity for larger k (smaller structures) is 
greater than the smaller k (larger structures). The modification of GR in the framework of 
f(R)- gravity, gives rise to an effective gravitational constant, G e ff, which is time and scale 
dependent parameter in contrast to the Newtonian gravitational constant. 

(vi) The linear density contrast relative to its value in a pure matter model, g(a ) = 5 m /a, 
for the all models starts from an early matter-dominated phase, g(a) = 1 , and decreases 
during history of the universe. 

(vii) The evolutionary behavior of the growth factor of linear matter density perturba¬ 
tions, f(z), shows that for the all models the growth factor for smaller k (larger structures) 
like the ACDM model fit the data very well. 


[1] M. Kowalski, et ah, (Supernova Cosmology Project), Astrophys. J. 686, 749 (2008). 

[2] E. Komatsu, et al., (WMAP Collaboration), Astrophys. J. Suppl. 192, 18 (2011); 

G. Hinshaw, et al., (WMAP Collaboration), Astrophys. J. Suppl. 208, 19 (2013). 

[3] P.A.R. Ade, et al., (Planck Collaboration), Astron. Astrophys. 571, A16 (2014). 

[4] M. Tegmark, et al., (SDSS collaboration), Phys. Rev. D 69, 103501 (2004); 

D.J. Eisenstein, et al., (SDSS collaboration), Astrophys. J. 633, 560 (2005); 

H. Lampeitl, et al., (SDSS collaboration), Mon. Not. Roy. Astron. Soc. 401, 2331 (2009). 

[5] T. Padmanabhan, Phys. Rep. 380, 235 (2003); 

P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75, 559 (2003); 

C.G. Tsagas, A. Challinor, R. Maartens, Phys. Rept. 465, 61 (2008); 


14 



M. Li, et al., Commun. Theor. Phys. 56, 525 (2011). 

[6] T.P. Sotiriou, V. Faraoni, Rev. Mod. Phys. 82, 451 (2010); 

A. De Felice, S. Tsujikawa, Living Rev. Relativ. 13, 3 (2010); 

S. Nojiri, S.D. Odintsov, Phys. Rept. 505, 59 (2011). 

[7] S. Nojiri, S.D. Odintsov, Phys. Rev. D 68, 123512 (2003). 

[8] S.A. Appleby, R.A. Battye, A.A. Starobinsky, JCAP 06, 005 (2010). 

[9] H. Motohashi, A. Nishizawa, Phys. Rev. D 86, 083514 (2012); 

A. Nishizawa, H. Motohashi, Phys. Rev. D 89, 063541 (2014). 

[10] Y. Sobouti, Astron. Astrophys. 464, 921 (2007). 

[11] K. Kararni, M.S. Khaledian, JHEP 03, 086 (2011); 

K. Karami, M.S. Khaledian, Int. J. Mod. Phys. D 21, 1250083 (2012). 

[12] A. Raccanelli, et al., Mon. Not. Roy. Astron. Soc. 436, 89 (2013). 

[13] D. Huterer, et al., Astropart. Phys. 63, 23 (2015). 

[14] M. Kunz, D. Sapone, Phys. Rev. Lett. 98, 121301 (2007). 

[15] E. Bertschinger, P. Zukin, Phys. Rev. D 78, 024015 (2008). 

[16] S. Baghram, S. Rahvar, JCAP 12, 008 (2010). 

[17] N. Mirzatuny, S. Khosravi, S. Baghram, H. Moshafi, JCAP 01, 019 (2014). 

[18] Y.S. Song, W. Hu, I. Sawicki, Phys. Rev. D 75, 044004 (2007). 

[19] T. Jacobson, Phys. Rev. Lett. 75, 1260 (1995). 

[20] R.G. Cai, S.P. Kim, JHEP 02, 050 (2005). 

[21] M. Akbar, R.G. Cai, Phys. Lett. B 635, 7 (2006); 

M. Akbar, R.G. Cai, Phys. Lett. B 648, 243 (2007). 

[22] K. Karami, M.S. Khaledian, N. Abdollahi, Europhys. Lett. 98, 30010 (2012). 

[23] G. Izquierdo, D. Pavon, Phys. Lett. B 639, 1 (2006). 

[24] H. Mohseni Sadjadi, Phys. Rev. D 73, 063525 (2006); 

H. Mohseni Sadjadi, Phys. Rev. D 76, 104024 (2007); 

H. Mohseni Sadjadi, Phys. Lett. B 645, 108 (2007). 

[25] J. Zhou, et al., Phys. Lett. B 652, 86 (2007). 

[26] Y. Gong, B. Wang, A. Wang, Phys. Rev. D 75, 123516 (2007); 

Y. Gong, B. Wang, A. Wang, JCAP 01, 024 (2007). 

[27] M. Jamil, M. Akbar, Gen. Relativ. Gravit. 43, 1061 (2011). 


15 



[28] A. Sheykhi, B. Wang, Phys. Lett. B 678, 434 (2009); 

A. Sheykhi, Phys. Rev. D 81, 104011 (2010). 

[29] N. Banerjee, D. Pavon, Phys. Lett. B 647, 447 (2007). 

[30] K. Karami, JCAP 01, 015 (2010); 

K. Karami, S. Ghaffari, Phys. Lett. B 685, 115 (2010); 

K. Karami, S. Ghaffari, Phys. Lett. B 688, 125 (2010); 

K. Karami, S. Ghaffari, M.M. Soltanzadeh, Class. Quantum Grav. 27, 205021 (2010). 

[31] K. Karami, et ah, JHEP 08, 150 (2011); 

K. Karami, A. Abdolmaleki, JCAP 04, 007 (2012); 

K. Karami, et al., Eur. Phys. J. C 73, 2565 (2013); 

K. Karami, et ah, Phys. Rev. D 88, 084034 (2013); 

A. Abdolmaleki, T. Najafi, K. Karami, Phys. Rev. D 89, 104041 (2014). 

[32] N. Radicella, D. Pavon, Phys. Lett. B 691, 121 (2010). 

[33] K. Barnba, C.Q. Geng, JCAP 11, 008 (2011). 

[34] A.A. Starobinsky, Phys. Lett. B 91, 99 (1980). 

[35] H. Motohashi, A.A. Starobinsky, J. Yokoyama, Prog. Theor. Phys. 123, 887 (2010); 

H. Motohashi, A.A. Starobinsky, J. Yokoyama, JCAP 06, 006 (2011). 

[36] S. Capozziello, et ah, Int. J. Mod. Phys. D 12 1969, (2003); 

S. Capozziello, V.F. Cardone, A. Troisi, Phys. Rev. D 71, 043503 (2005). 

[37] N.J. Poplawski, Phys. Lett. B 640, 135 (2006); 

S. Capozzielo, et ah, Astrophys. Space Sci. 342, 155 (2012); 

V.F. Cardone, S. Camera, A. Diaferio, JCAP 02, 030 (2012). 

[38] J.c. Hwang, H. Noh, Phys. Rev. D 71, 063536 (2005). 

[39] S. Tsujikawa, Phys. Rev. D 77, 023507 (2008). 

[40] A.A. Starobinsky, JETP Lett. 86, 157 (2007). 

[41] S. Tsujikawa, Phys. Rev. D 76, 023514 (2007); 

S. Tsujikawa, K. Uddin, R. Tavakol, Phys. Rev. D 77, 043007 (2008). 

[42] H. Motohashi, A.A. Starobinsky, J. Yokoyama, Int. J. Mod. Phys. D 18, 1731 (2009). 

[43] P.J.E. Peebles, Principles of Physical Cosmology , Princeton University Press (1993); 

T. Padmanabhan, Structure formation in the universe , Cambridge University Press (1993). 

[44] J.D. Bekenstein, Phys. Rev. D 9, 3292 (1974). 


16 



[45] R.M. Wald, Phys. Rev. D 48, 3427 (1993); 

G. Cognola, et al., JCAP 02, 010 (2005). 

[46] L.G. Jaime, L. Patino, M. Salgado, arXiv: 1206.1642. 

[47] W. Hu, I. Sawicki, Phys. Rev. D 76, 064004 (2007). 

[48] K. Bamba, C.Q. Geng, C.C. Lee, JCAP 08, 021 (2010). 

[49] A.D. Dolgav, M. Kawasaki, Phys. Lett. B 573, 1 (2003). 

[50] V. Faraoni, Phys. Rev. D 74, 104017 (2006). 

[51] S. Tsujikawa, et al., Phys. Rev. D 80, 084044 (2009). 

[52] A. de la Cruz-Dombriz, A. Dobado, A.L. Maroto, Phys. Rev. Lett. 103, 179001 (2009). 

[53] H. Motohashi, A.A. Starobinsky, J. Yokoyama, Int. J. Mod. Phys. D 20, 1347 (2011). 

[54] K. Bamba, C.Q. Geng, C.C. Lee, JCAP 11, 001 (2010); 

K. Bamba, C.Q. Geng, C.C. Lee, Int. J. Mod. Phys. D 20, 1339 (2011). 

[55] S.A. Appleby, R.A. Battye, Phys. Lett. B 654, 7 (2007); 

S.A. Appleby, R.A. Battye, JCAP 05, 019 (2008). 

[56] S. Capozziello, et al., Phys. Rev. D 84, 043527 (2011). 

[57] L. Verde, et al., Mon. Not. Roy. Astron. Soc. 335, 432 (2002); 

E. Hawkins, et al., Mon. Not. Roy. Astron. Soc. 346, 78 (2003); 

E.V. Linder, Astropart. Phys. 29, 336 (2008). 

[58] C. Blake, et al., Mon. Not. Roy. Astron. Soc. 415, 2876 (2011). 

[59] R. Reyes, et al., Nature 464, 256 (2010). 

[60] M. Tegmark, et al., Phys. Rev. D 74, 123507 (2006). 

[61] N.P. Ross, et al., Mon. Not. Roy. Astron. Soc. 381, 573 (2007). 

[62] L. Guzzo, et al., Nature 451, 541 (2008). 

[63] J. da Angela, et al., Mon. Not. Roy. Astron. Soc. 383, 565 (2008). 

[64] P. McDonald, et al., Astrophys. J. 635, 761 (2005). 


17 













































































Starobinsky f(R) Model Starobinsky f(R) Model 



Starobinsky f(R) Model Starobinsky f(R) Model 



FIG. 6: The variations of RFr, the screened mass function G e s/G, the linear density contrast 
relative to its value in a pure matter model g = 5/a and the growth factor f(z), versus redshift z 
for the Starobinsky model. 
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FIG. 7: Same as Fig. 6 but for the Hu-Sawicki model. 
















FIG. 8: Same as Fig. 6 but for the Exponential model. 















FIG. 9: Same as Fig. 6 but for the Tsujikawa model. 















FIG. 10: Same as Fig. 6 but for the AB model. 














TABLE I: The observational data for the linear growth rate f 0 bs(z)- 

z 0.15 0.22 0.32 0.35 0.41 0.55 0.60 0.77 0.78 1.4 3.0 

/obs 0.51 0.60 0.654 0.70 0.70 0.75 0.73 0.91 0.70 0.90 1.46 

la 0.11 0.10 0.18 0.18 0.07 0.18 0.07 0.36 0.08 0.24 0.29 
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